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Using the spin-polarized tunneling model and tak- 
ing into account the basic physics of ferromagnetic semi- 
conductors, we study the temperature dependence of the 
tunneling magnetoresistance (TMR) in the diluted mag- 
netic semiconductor (DMS) trilayer heterostructure system 
(Ga,Mn)As/AlAs/(Ga,Mn)As. The experimentally observed 
TMR ratio is in reasonable agreement with our result based 
on the typical material parameters. It is also shown that the 
TMR ratio has a strong dependence on both the itinerant- 
carrier density and the magnetic ion density in the DMS elec- 
trodes. This can provide a potential way to achieve larger 
TMR ratio by optimally adjusting the material parameters. 



PACS numbers: 72.25.-b, 73.40.Ty, 73.40.Gk, 75.50.Pp 



I. INTRODUCTION 

Since the early experiment done by Jullierejl] there 
has been considerable interest in studying the tunnel- 
ing magnetoresistance (TMR) effect in magnetic tunnel 
junctionsou due to its promising applications in digital 
storage and magnetic sensor technology.^ The conven- 
tional electrodes of magnetic tunnel junctions are made 
up of the ferromagnetic transition metals such as Fe, Co, 
Ni, and their alloys. Usually the TMR effect in the con- 
ventional magnetic tunnel junctions is described by the 
Julliere modelB'B The electrons near at the Fermi levels 
of the electrodes participate in the tunneling process and 
have the main contribution to the tunneling conductance 
and so to the TMR effect. 

Recent discovery of the ferromagnetism in the di- 
luted magnetic semiconductors (DMS) has opened the 
rich field of an interesting intemiay between magnetic 
and semiconducting properties ,En3 which makes it pos- 
sible to combine the information processing and data 
storage in one material. At high concentration c of 
randomly distributed Mn 2+ ions doped in GaAs sam- 
ples with high hole density c*, Gai_ x Mn x As compounds 
exhibit ferromagnetism with transition temperature as 
high as 110 K at certain value of x. The origin of 
ferromagnetism in the DMS can be explained by the 
itinerant holes through the p-d exchange coupling be- 
tween the valanc&iband holes and the spins of im- 
purity Mn ionsJa £l In order to investigate the spin- 
polarized electron transport phenomenon in the DMS fer- 
romagnetic semiconductor heterostructures, the epitax- 
ial tunnel junctions have been successfully fabricated in 



the form of Ga .96iMno.o39As/AlAs/Gao.96iMn .o39As. 
The TMR effect was observed and the maximum 
TMR xatio thus obtained is as high, as 44.4% at 
4.2 K.Q The subsequent experimentt3 also realized 
the TMR effect in the DMS tunnel junctions of 
Gao.95Mno.05As/AlyGa1_yAs/Gao.97Mno.03As. 

The-,,-, |— 1 |— J Julliere 
model- - with several necessary extensions-irE-l is known 
to describe the tunnel junctions, in which the two elec- 
trodes are made of the conventional ferromagnetic metals 
or the half-metallic manganites such as Lai-^Sr^MnOs. 
For the manganite tunnel junctions, various tunneling 
processes are included in the model as follows. The elas- 
tic spin-flip scattering of the spin polarized carriers by 
the local magnetic impurities within the insulating bar- 
rier and the inelastic spin-flip scattering induced by the 
thermal collective excitations of local spins at the inter- 
faces between the insulating barrier and the manganite 
electrodes are the two basic ingredients. By taking into 
account the temperature dependence of the spin polar- 
ization of itinerant carriers, one can calculate the tem- 
perature dependence of the TMR ratio.Ej In the paper, 
we have applied the above generalized Julliere model to 
the DMS tunnel junctions and obtained the various prop- 
erties of the TMR effects of the system. 

The paper is organized as follows. In Sec. II, we pro- 
pose a spin-polarized tunneling model for the DMS tun- 
nel junctions by generalizing the previous spin-polarized 
tunneling model in order to correctly take into account 
the basic physics of the DMS. In Sec. Ill, we show that 
the experimentally observed TMR ratio in the DMS tun- 
nel junctions can be well reproduced from our calcula- 
tions with the appropriate choice of typical material pa- 
rameters and also show that the TMR ratio thus obtained 
is strongly dependent on both the itinerant-carrier den- 
sity and the magnetic ion density in the DMS electrodes. 
Finally we conclude with a brief summary in Sec. IV. 



II. SPIN-POLARIZED TUNNELING MODEL 

We will briefly summarize the basic physics of ferro- 
magnetism occurring in certain DMSLcompounds follow- 
ing the seminal work of Konig et al.n 

In the doping range of our interest, that is, 0.035 < 
x < 0.06, Gai^Mn-rAs can be described as a degenerate 
semiconductor with relatively high hole concentrations 
and the Mn 2+ ions act as local magnetic impurities. 

The Hamiltonian describing the basic physics of the 
DMS can be written as 
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H= /"d3 r ^^t (r)( _^_v 2 -M*)^(r) 

J a 

+ J pd Jd 3 rS(r)- S (r), (1) 

where S(r) stands for the spin density of the local mag- 
netic impurities, which is given by J2i S/<5(r — Rf ) with 
R/ denoting the positions of magnetic impurities. The 
field operator ^(r) (^/ cr (r)) is the creation (annihila- 
tion) operator for the itinerant hole at position r with 
spin <7. The s(r) represents the spin density of itinerant- 
carriers, which can be written in terms of field operators 
as ^X)aa' *3r( r ) 1 CTcr' v I > cr'(r), where t are the three Pauli 
spin matrices. The m* and fi* are the effective mass 
and chemical potential of itinerant carriers, respectively. 
The J p d stands for the p-d exchange coupling strength 
between itinerant holes and Mn 2+ ion impurity spins. 
The recent experiments suggest that the itinerant hole 
density c* is much lower than the magnetic ion density 
c satisfying c*/c <C 1. Here we focus our attention on 
the DMS with ferromagnetic ground state. Using the 
Holstein-Primakoff transformations and performing the 
coarse graining procedures, the local spin density S(r) 
can be written in terms of the newly defined bosonic fields 
b'(r),b(r) as follows 

5+(r) = yj2cS-tf(r)b(r)b(r), (2) 

S- (r) = fet ( r ) ^2c5-6t(r)6(r), (3) 

S z (r) = cS-tf (r)6(r). (4) 

The self-consistent spin wave approximation developed 
in Ref. || yields the following effective Hamiltonian for 
the DMS system 

Heff = ^ EkffC^Ckff + ^ ^qbq&q, (5) 

kcr |q|<<3c 

where Eko- stands for the band energy of the itinerant 
holes with spin a — ±1, which is given by 

e-k a = £k + y a ( 6 ) 

and £k = h 2 k 2 /(2m*). The spin wave dispersion f2 q is 
given by 

fi q = ^f-( n l- n \) - H [/( £ kl) - /(£k+qt) 

k 



where Ck and 6 q are the itinerant hole and the spin wave 
annihilation operators in the momentum space respec- 
tively, and f(x) is the Fermi-Dirac distribution func- 
tion. Here nj_ (nj) is the spin-down (spin- up) itinerant 
hole density satisfying nj + n± — c* , the exchange gap 
A = cJpdS, and q c is the Debye cutoff for the spin waves 
of local impurity spins with q 3 = 6tt 2 c. At finite tem- 
peratures, the spin wave dispersion can be generalized 
by imposing the following self-consistency conditions for 
the finite temperature exchange gap : A(T) = J p d(S z ). 
Here (S z ) represents the thermal average of the impurity 
Mn ion spins, which are approximately calculated by the 
following formula 

q|<<?c 

= y {S-n B (n ci ) + (2S+l)n B [(2S+m (1 }}, 

q|<9c 

(8) 

where B s (x) is the Brillouin function, hb(i) the Bose- 
Einstein distribution function, and [3 = 1/ksT. The sec- 
ond term in the second equality of Eq. (8) describes 
how the thermally induced spin waves from the inde- 
pendent Bose statistics with no limit in the number of 
spin waves indeed reduce the magnetization of the sys- 
tem. The third term takes into account the correct spin 
kinematics, which rules out the unphysical states. 

Based on the above scheme, one can calculate the mag- 
netization of the Mn impurity spins and the band struc- 
ture of the itinerant holes approximately self-consistently 
by solving a set of coupled equations. In this theoretical 
framework, the non-monotonic dependence of the Curie 
temperature on the itinerant hole density c* has been 
successfully obtained.^ 

Now we consider the DMS magnetic tunnel junctions 
in which both electrodes are made of the ferromagnetic 
DMS Gai-^Mn^As compounds. The carrier tunneling 
through a thin AlAs insulating barrier involves removing 
a carrier at one side and creating it at the other side with 
spin conserving or spin-flipping. The spin-flip tunneling 
processes can be divided into two different categories: 
The elastic spin-flip tunneling arising from the Mn 2+ ion 
impurity spins located within the insulating barrier and 
the inelastic spin-flip tunneling induced by the thermal 
collective excitations of the impurity spins near at the 
interfaces between the insulating barrier and the DMS 
electrodes. 

For the collinear magnetization configurations, the 
tunneling Hamiltonian of the DMS tunnel junctions can 
be written as 

H = H L + H R + H T , (9) 

# L = £>p- c p- c p- + Yl Q ^IlK,l, (10) 
p°" |q|<<?c 
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J kp a kcr L pc r 



|q|<» 



(11) 



kp a kcr C P--cr 



H.c. 



|T"| 2 , respectively. By making a transformation for the 
tunneling Hamiltonian from the collinear magnetization 
configurations to arbitrary relative angle 9 between the 
spin orientations of the two electrodesEZI and using the 
standard Green's function methods Jla one can obtain the 
tunneling conductance G(9) at zero bias up to the leading 
order in the tunneling amplitudes, which is given by 



E { r kpq[ 5 i(q) + ^(q)] (4t c pt - 4^) 



kp,|q|<9c 



G(9) = — \T°\ 2 [ Pi (e F ) + PT (e F )] 2 [l + P 2 cos9 



-H.c.} + -^= Yl { T k Pq [v / 2^^L + 



kp.|q|<(j c 



(4| c pT 



C P l d kTj 



H. 



} 



(12) 



Here Hi (Hn) is the effective Hamiltonian for the bulk 
DMS in the left (right) electrode, which has been derived 
from the original Hamiltonian of Eq. (1) as explained 
earlier. The operator dko- (c po -) is the carrier annihilation 
operator with momentum k (p) and spin a of the right 
(left) DMS electrode, and S*(q) = cS a - & qa 6 qQ with 
a = L, R. 

The term with coefficient T k in Hi represents the 
spin-flip tunneling arising from the Mn 2+ ion impurity 
spins inside the AlAs insulating barrier. This type of 
spin-flip tunneling has been introduced to the tunneling 
Hamiltonian for the conventional ferromagnetic, jtunncl 
junctionst 2 ] and the manganite tunnel junctions.c3'c3 The 
imaginary i for the spin-flip tunneling matrix elements in- 
dicates that the elastic spin-flip tunneling processes are 
incoherent in any magnetization configurations. As a 
generalization of the typical tunneling HamiltoniantLj to 
the DMS tunnel junctions, the last two terms have been 
included in the tunneling Hamiltonian of Eq. (12), which 
describe the inelastic spin-flip tunneling process induced 
by the thermal spin wave collective excitations of local 
spins near at the interfaces between the insulating barrier 
and the DMS electrodes. Since this kind of carrier tunnel- 
ing involves the emission or absorption of the local spin 
collective excitations near the left or right interfaces, it is 
an inelastic tunneling process. The inelastic spin-flip tun- 
neling may occur not only directly at the interfaces but 
also within the coherence length of the interfaces, which 
is relatively long comparing to the barrier width. It then 
allows us to use the bulk spin wave dispersion relation 
in the tunneling Hamiltonian. The fact that the wave 
vector k is independent of p in Eq. (12) implies that we 
only consider the incoherent tunneling processes, which 
are known to be appropriate for the defect-populated 
insulating barriers and the rough interfaces. As usual, 
we make the following assumptions that the strengths 
of various tunneling amplitudes, that is, |T^ p | 2 , |T£ p | 2 , 
and 1 7k' | 2 are independent of momenta k, p, q and then 
will be taken to be their average values |T°| 2 , |T'| 2 , and 



+ 7 (1 -P 2 cos0)], 



(13) 



where P stands for the itinerant carrier spin polarization 
given by 



P 



pi(sf) - gt(gg) 
PI&f) +p^{e F y 



(14) 



and Pi(sf) (pt( £ ^)) i s the density of states (DOS) of 
the spin-down (spin-up) carriers at the Fermi level of the 
DMS electrodes. The relative angle 9 may be tunable by 
applying a magnetic field and 7 is given by the following 
relation 



1 = ll +Ar ] cS{cS~(S z )), 



(15) 



2 and 7] 

| T 0|2 



T"| 2 /|T| 2 with |T| 2 = 
2| T"| 2 is usxujilly 



since c 



where 71 = \T'\ 2 /\T 
\T°\ 2 + 2c 2 S 2 \T"\ 2 , , 
about two orders of magnitude smaller than |r°| 2 .t3 In 
deriving Eq. (15), we have used the following relation 
V- 1 E| q |< gc ( fo q & q> = cS ~ ( sz ) ■ The corresponding tun- 
neling resistance is given by R(9) — 1/G(6). 

The TMR ratio is defined by the following formula 



AR 



Rap — Rp 
Rp~ ' 



(16) 



where Rap and Rp represent the tunneling resistance in 
the antiparallel and parallel spin alignments of the two 
ferromagnetic electrodes respectively, and the maximum 
TMR ratio for the DMS tunnel junctions can be obtained 



AR 
~Rp~ 



2(1 -l)P 2 



1-P 2 + 7(1 + P 2 )' 



(17) 



If the electrodes of the DMS tunnel junctions are 
made of the non-degenerate DMSs, the tunneling conduc- 
tance will be proportional to the itinerant carrier density 
and the TMR ratio will be more appropriately given by 
Eq. (17) with the carrier spin polarization replaced by 
P = (ni — n^)/(ni + n-p). The different definitions of 
the carrier spin polarization in the TMR ratio for the de- 
generate and non-degenerate DMS tunnel junctions are 
due to the distinctly different transport properties in the 
corresponding DMSs. 
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III. RESULTS AND DISCUSSIONS 

In Fig. 1 and Fig. 2 are illustrated the temperature 
dependences of the impurity Mn ion magnetization of 
the DMS, the itinerant-carrier spin polarization, and the 
maximum TMR ratio for the DMS tunnel junctions at 
several different itinerant-carrier densities c* with the 
following choice of material parameters: 71 = 0.692, 
77 = 0.005 nm 6 , and m* = 0.5m e , J pd = 0.15 cY.nm 3 
appropriate for the bulk Gai-^Mn^As compoundsEl and 
fixed impurity Mn ion density c = 1 nm -3 , where m e is 
the free electron mass. Figure 1 (a) shows the impurity 
Mn ion magnetization (S z ) normalized by c. One can no- 
tice that the critical temperature T c decreases with the 
further increase of the itinerant carrier density c* from 
c* =0.1 nm~ 3 . From the solid curves of Fig. 1 (b) cor- 
responding to c* =0.1 nm~ 3 , one can see the following 
features that at relatively low temperatures, the spin-up 
carrier band is completely lifted from the spin-down car- 
rier band by the large exchange gap A(T). At the mean- 
field level, this leads to the complete spin polarization of 
the itinerant carriers, i.e. P = 100%, exhibiting the half- 
metallic feature. In this case the TMR ratio is reduced to 
(AR/R p ) max = I/7— 1, which decreases with the increase 
in the population of the thermally excited spin wave col- 
lective excitations leading to the increase of 7 as shown in 
Fig. 1 (c). As the temperature increases, the size of the 
exchange gap A(T) is gradually reduced. With the fur- 
ther increase of temperature beyond a certain value, the 
number of the majority spin-down carriers begins to de- 
crease, while that of the minority spin-up carriers starts 
to increase. At high temperatures, the carrier spin polar- 
ization decreases rapidly with temperatures, which makes 
the TMR effect disappear above a critical temperature. 
The factor of (1 - 7) in Eq. (17) indicates that the TMR 
effect may vanish below the critical temperature, which is 
due to the competition between the spin-conserving tun- 
neling amplitude |T°| 2 and the spin- flip tunneling ones 

\t'\ 2 ,\t"\ 2 . 

For higher carrier densities, say, c* = 0.2 nm -3 , the 
Fermi energy of the system increases and then both spin 
states are occupied losing the half-metallic feature as 
shown in Fig. 1 (b). Since the partial carrier spin polar- 
ization is not beneficial to the TMR effect, it leads to the 
low TMR ratio. Figure 2 (a) illustrates the impurity Mn 
ion magnetization (S z ) normalized by c at different val- 
ues of c*. In contrary to Fig. 1 (a), one can see that T c de- 
creases with the decrease of the itinerant carrier density 
c* from c* = 0.05 nm -3 exhibiting the non-monotonic 
dependence of T c on c*, which has been emphasized in 
Ref. ^. For c = 1 run -3 , the optimal carrier densities, 
which will achieve the maximum TMR ratio, are shown 
to be in the range of 0.05 nm -3 < c* < 0.1 nm~ 3 . 

Recently Hayashi et al. have performed an experi- 
ment for Gao.96iMno.o39As/AlAs/Gao.96iMno.o3gAs tun- 
nel junctions and obtained the TMR ratio of about 44% 
at 4.2 Kfl By fitting the above experimental data to 



our calculations, we were able to estimate the approx- 
imate values of material parameters: c* = 0.13 nm -3 , 
c = 0.86 nm~ 3 , and 71 = 0.692. The above value of car- 
rier density c* is about 15 % of the Mn ion concentration 
c. The value of 71 indicates that there exists strong spin- 
flip scattering by the Mn ions within the insulating bar- 
rier. During the fabrication process of the DMS tunnel 
junctions, the highly doped Mn ions in the Gai-^Mn^As 
layers can easily diffuse into the AlAs insulating layer. 
Reducing the Mn ions in the insulating layer is crucial to 
achieve a large TMR ratio in the DMS tunnel junctions. 

Figure 3 shows the temperature dependences of the im- 
purity Mn ion magnetization of the DMS, the itinerant- 
carrier spin polarization, and the maximum TMR ra- 
tio for the DMS tunnel junctions at several values of 
impurity Mn ion densities c with fixed carrier density 
c* =0.1 nm~ 3 . Figures 3 (a) and (b) demonstrate that 
the carrier band structure is mainly determined from the 
size of the exchange gap A(T). At low impurity densi- 
ties, say, c = 0.6 nm~ 3 , the strength of the ferromagnetic 
coupling is weak, which induces partial carrier spin po- 
larization. By increasing c, the ferromagnetic coupling 
becomes strong, which completely polarizes the carrier 
spins above a certain value of c. So achieving the high 
doping of the impurity Mn ions is necessary to have much 
higher T c and larger TMR ratio at high temperatures. 

When the magnetic ion density is beyond a certain 
large value, say, c = 1.5 nm~ 3 , the impurity ion magne- 
tization and the itinerant-carrier spin polarization seem 
to exhibit the behaviors of the first order ferromagnetic 
transition as shown in Fig. 3 (a) and (b). In this case, 
the first order ferromagnetic transition is induced by the 
p-d exchange coupling, in contrast to that for the two- 
dimensional quantum well system at very low carrier den- 
sities, which is acgued to be due to the carrier-carrier 
correlation effectsJlS 

Our results clearly demonstrate that the TMR effects 
in the DMS tunnel junctions are dominated by both 
the ferromagnetic coupling strength and the carrier band 
structure through the impurity ion density and the itin- 
erant carrier density in the DMS electrodes. The elastic 
and inelastic spin-flip tunnelings play an important role 
in reducing the TMR ratio. 

IV. SUMMARY 

In summary, we have presented the TMR formula for 
the DMS tunnel junctions by generalizing the previous 
spin-polarized tunneling model and taking into account 
the basic physics of the DMS. The experimentally ob- 
served TMR ratio was reproduced by the appropriate 
choice of the typical material parameters. It was also 
shown that the TMR ratio has a strong dependence on 
the itinerant-carrier density and the magnetic ion den- 
sity in the DMS electrodes as well. In view of the re- 
cent predictions of the room temperature ferromagnetic 
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semiconductors .0 the DMS tunnel junctions can be prac- 
tically very useful in the future digital storage and mag- 
netic sensor technology. Our analysis of the TMR ratio 
can be potentially useful to achieve larger TMR ratio by 
optimally adjusting the material parameters. 
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Fig. 1. (a) Impurity Mn ion magnetization of the DMS, 
(b) itinerant-carrier spin polarization P, and (c) the max- 
imum TMR ratio in the DMS tunnel junctions as a func- 
tion of temperature for several different itinerant-carrier 
densities c* = 0.3, 0.2, 0.18, 0.15, 0.1 nm" 3 with the 
following choice of material parameters m* — 0.5m e , 
J pd = 0.15 eV nm 3 , 71 = 0.692 and 77 = 0.005 nm 6 , 
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and the fixed value of c = 1 nm 3 . 



and the fixed value of c = 1 nm 3 . 
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Fig. 2. (a) Impurity Mn ion magnetization of the DMS, 
(b) itinerant-carrier spin polarization P, and (c) the max- 
imum TMR ratio in the DMS tunnel junctions as a func- 
tion of temperature at several different itinerant-carrier 
densities c* = 0.05, 0.03, 0.01, 0.001 nm" 3 with the 
following choice of material parameters m* — 0.5m e , 
J pd = 0.15 eV nm 3 , 71 = 0.692 and r\ = 0.005 nm 6 , 



Fig. 3. (a) Impurity Mn ion magnetization of the DMS, 
(b) itinerant-carrier spin polarization P, and (c) the max- 
imum TMR ratio in the DMS tunnel junctions as a func- 
tion of temperature for several different impurity Mn 
ion densities c = 0.6, 0.7, 0.9, 1.2, 1.5 nm~ 3 with the 
following choice of material parameters m* = 0.5m e , 
J pd = 0.15 eV nm 3 , 71 = 0.692 and r\ = 0.005 nm 6 , 
and the fixed value of c* =0.1 nm -3 . 
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